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' We report on the stationary dynamics in classical Sinai billiard (SB) 



corresponding to the unit cell of the periodic Lorentz gas (LG) formed by 



' square lattice of length L and dispersing circles of radius R placed in the cen- 



ter of unit cell. Dynamic correlation effects for classical particles, initially 
^ ' distributed by random way, are considered within the scope of deterministic 

and stochastic descriptions. A temporal analysis of elastic reflections from 
the SB square walls and circle obstacles is given for distinct geometries in 
terms of the wall-collision and the circle-collision distributions. Late-time 
steady dynamic regimes are explicit in the diffusion exponent z{R), which 
plays a role of the order-disorder crossover dynamical parameter. The bal- 
listic {zq = 1) ordered motion in the square lattice {R = 0) switches to 
the superdiffusion regime with zi = 1.5, which is geometry-independent 
when R < L\/^/4:. This observed universal dynamics is shown to arise from 
long-distance particle jumps along the diagonal and nondiagonal Bleher cor- 



1 



ridors in the LG with the infinite horizon geometry. In the corresponding 
SB, this universal regime is caused by the long-time wall-collision memory 
effects attributed to the bouncing-ball orbits. The crossover nonuniversal 
behavior with 1.5 < z < 2 is due to geometry with L\/2/4 < R < L/2, 
when only the nondiagonal corridors remain open. All the free-motion 
corridors are closed in LG with finite horizon (R > L/2) and the inter- 
play between square and circle geometries results in the chaotic dynamics 
ensured by the normal Brownian diffusion {z2 = 2) and by the normal 
Gaussian distribution of collisions. PACs: 05.45.Gg , 05.40.Fb, 45.50.Tn. 
Key words: Sinai billiard, Lorentz gas, collision statistics, anomalous dif- 
fusion, order-disorder crossover, chaos. Corresponding author Valery B. 
Kokshenev, valery@fisica.ufmg.br 
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I. INTRODUCTION 



A mathematical exploration of two-dimensional Lorentz gas (LG) with periodic config- 
uration of scatterers, related to dispersed billiards, was initiated [1] by Bunimovich and 
Sinai. The process of establishing equilibrium dynamics for randomly distributed, noninter- 
acting classical particles driven by elastic collisions with the billiard walls and the scatterer 
of fixed geometry are commonly discussed in terms of the particle velocity- velocity and/or 
displacement-displacement correlation functions, which are related to a diffusion coefficient 
through the Einstein-Green- Kubo formula (see e.g. Ref. [2]). Besides the ergodicity, the 
entropy, the Lyapunov exponent, the mixing property, and among other interesting physical 
observables in LG with finite and infinite horizons, the enhanced diffusion studied numeri- 
cally through the collision distribution function, has been expressed as of great importance in 
Ref. [3] . This stimulated our subsequent investigations of the coUision statistics for nonescape 
particles observed through their survival probability in weakly open chaotic (Sinai billiard [4] , 
(SB)) [5] and non-chaotic (circle and square biUiards) [6] classical systems. Remarkably, col- 
lision statistics provided new insight into the delicate mathematical problem of the interplay 
between regular and irregular segments of the billiard boundary, as has been demonstrated 
for the case of the so-called almost-integrable systems, presented by open [7] and closed [8] 
rational polygons. In the current study we focus on the stationary superdiffusion dynamics 
in the closed SB. Analysis is given through the billiard- wall and the scatterer-disk random 
collision statistics and is based on the dynamical correspondence between SB and LG. Some 
of the findings of this study were preliminary communicated in Ref. [8] . 

The dynamics of particles (of unit mass and unit velocity) moving in two-dimensional 
closed region (billiard table) and dispersed by obstacles is governed by billiard-boundary 
geometry. The SB with the square-wall table of length L and the disk with radius R can 
be formally treated as the unit cell of a periodic LG, which is the two-dimensional periodic 
crystal formed by a regular set of circular scatterers (of radius R) centered at distances L. 
This implies that SB and LG are dynamically equivalent classical systems. This statement 
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is justified by a step-to-step correspondence that can be established between a certain orbit 
in the closed SB and the corresponding trajectory in LG. For the case of the low-density LG, 
given by geometry with R < L/2, this correspondence is exemplified in inset A in Fig. 1. As 
seen for this LG configuration, there exist trajectories in which particles never collide with 
scatterers. Those particles, which move through such unbounded trajectories, have therefore 
an infinite horizon [2,3]. In turn, the free-motion trajectories belong to infinite corridors [2]. 
As shown by Bleher [2], the principal corridors are open for small disk radii limited by 
< i? < L-\/2/4. This is illustrated by the diagonal and nondiagonal corridors in inset B 
in Fig.l. When the scatterer radius achieves the magnitude of L\/2/4, the last nondiagonal 
corridor closes. With further increasing the radius, the principal corridors disappear at L/2 
and scatterers start to overlap lattice cells. The particle trajectories become bounded thus 
having a finite horizon (see inset C in Fig.l). 

In the infinite-horizon LG, the enhanced-diffusion motion regime was theoretically pre- 
dicted in Ref . [2] , through the asymptotic statistical behavior of the late-time particle dis- 
placements. But no description for the evolution of motion regimes with geometry was given 
(see also Ref. [9]). We therefore reformulate the problem of random particle displacements 
into that of random collisions. This yields a description of the late-time stationary dynamics 
in LG, as well as in the corresponding SB, through the diffusion exponent z{R), which is 
continuous with R. The paper is organized as follows. In Sec. II we develop the billiard 
collision statistics based on the alternative deterministic and stochastic approaches. Also, 
we introduce the collision distribution function given in terms of the dynamic observables, 
which are available in simulation of SB with finite and infinite horizon geometries in the 
corresponding LG. Conclusion is drawn in Sec. III. 
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II. COLLISION STATISTICS 



A. Distribution Function and Dynamic Characteristics 

The trajectories of classical particles of unit mass moving with unit velocities in SB (or 
in LG) are preserved by the Liouville measure [10], namely 

(i/x(x) = ^ ^ dxdydO. (1) 

This is introduced in the phase space through the billiard table area A, given by the variable 
scatterer radius R and the fixed boundary side L. The coordinate set x = (x, y, 9) includes 
the particle position and the velocity launching angle 9 — [0,27r], which is counted of the 
X-axis of the square billiard table. The mean collision time Tc{R), which is due to the 
two consequent elastic random collisions with the boundary wall or with the scatterer, are 
respectively defined [10] by 

ri-Hi^) = ^orrW(i?) = ^. (2) 



Here Py, — 4L and Pg — 27rR are corresponding collision perimeters of the table with 
accessible area A = L'^ — txR^ ( for i? < L/2). 

The collision distribution function D{n, t) is a probability of a particle to collide n times 
with the fixed billiard boundary within a time t (for rigorous definition see e.g. Ref. [3]). 
Also, it can be introduced through the billiard mean collision-number equations, namely 

f f°° t 

Ucit) =< n(x,t) >c= / n(x,t) (i/i,(x) = / n D{n,t)dn = — , (3) 

J Jo 

where x stands for the boundary (wall and/or scatterer) position set and Tc{R) is given in 
Eq.(2). This results in 

(i/i[x(n, t)] 



D{n, t) 



(4) 



dn 

where the Liouville measure, defined in Eq. (1), is given by the inverse function to n(x,t). 
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The distribution (4) provides a rich information on the boundary-memory effects in 
chaotic SB with a fixed geometry. In other words, dynamic correlation effects can be charac- 
terized by the nonzero central moments of order m = 2, 4... defined for the random numbers 
n as 

/•oo 

A"^nc(i)=/ [n-nc{t)]"'D{n,t)dn. (5) 
Jo 

A"^nc{t) describes m-order deviation from the mean collision number nc{t) given in Eq.(3). 
Basically, we focus on low-order dynamic correlation effects presented by the variance of the 
random collision numbers n 

AVW = ^(AV)e, (6) 

T 

expressed through the variance for particle displacements (A^r)c. This fundamental equation 
will be deduced below using the orbit-trajectory correspondence visualized in Fig. 1. Also, 
Eq. (6) permits one to employ the well known temporal variance- displacement equation (see 
e.g. Ref. [11]) 



(A^r)e , with t » Te. (7) 



In this way, the diffusion exponent z{R) introduces a description for distinct stationary 
regimes in billiards with different R. For particles of unit velocity, the mean free path 
ic — Tc can be specified for the wall and the scatterer coUisions with the help of Eq.(2). 

The wall-collision and the scatterer- collision statistics is given by the corresponding 
distribution functions D^^\n, t) and D^^\n, t). These are defined by the equations n^^\t) — 
t/rc^^ and n^c\t) — t/rc^\ which extend Eq.(3), as well as Eqs. (4) and (5). Furthermore, 
from the same equations one has ri^'^'ni"''' = rjf'^n^c^ = t that, with accounting for the total 
number of collisions nc{t) — n'^^ -\- n'f ^ with TcUc — t , results in the overall mean collision 
frequency, namely 

Here the mean wall (ri"''*) and scatterer (ri*'*) collision times are defined in Eq.(2). 
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B. Wall Collisions in Square Billiard 



If one ignores spliting effects caused by tt/ 2- angle vertices [7], particle in the closed square 
billiard are subjected to an ordered orbit motion driven by billiard walls. This imphes that 
the velocity launching angle 9 is the integral of motion for a given orbit, which can be 
therefore presented by a straight-line trajectory in the corresponding LG (see inset D in Fig. 
1). A whole number of intersections of this line with the unit-cell boundaries, encounted 
in X and y directions in a time corresponds to the following square-billiard, wall-collision 
number, namely [12] 

nj,"')(^,i) = l(cos^ + sin^), (9) 
Ij 

established for a given orbit by the velocity angle 9. In view of the point symmetry of 
square lattice, the angle domain is reduced to0<^<7r/4. Extending Eq.(3) for the 
wall-ordered motion, the corresponding characteristic time t^\9) follows from equation 

ni;"\9,t) = t/i^\9). This gives 

1 cos ^ -I- sin ^ . . 

WW) ^ ^ ' 

obtained with the help of Eq.(9). Consequently, the mean collision number is 

7r/4 A 

/ (w),^ ^\ 4i f cos9 + sm9 f , , 4i t 
n,o(t) = {nr{0,t))^ = — y d9 J dxdy = — = (11) 

Here the averaging procedure is elaborated over all equivalent angles 9. In turn, Eq.(ll) 
defines the mean collision time Tco that agrees with Eq.(2), where A — L'^ and Pyj — 4L. 
The variance for the random wall-collision number is 

AVo(i) = {[n^^\0,t) - n^{t)r)^^ = (l^ + i " ^'o, (12) 

with is obtained in Eq.(ll). The function 9{n, t) — 0.5 arcsin [{ATcon/'rrt)'^ — l] , inverse to 
the function n{9,t) given in Eq.(9), results in the collision distribution function DQ{n,t) — 
47r~^ \d9{n,t)/dn\ defined in Eq. (4). A straightforward estimation for the wall-collision 
distribution yields 
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.{w)i ,x _ 16 . -1 / TT 1 



^0 t) = —r, 7k 7 - o «cr sin 



'1^)2 _ 1 



for 7r/4 < n/rico < 7rV2/4, otherwise Bf^{n, t) = 0. (13) 

One can verify that Eqs. (11), (12) and (13) are selfconsistent hence obey Eqs.(3) and (5). 
As seen from Eqs. (7) and (12), the dynamical diffusion exponent zq — 1, i.e., the dynamic 
regime in the square bilhard is balhstic. Meantime, balhstic trajectories were employed to 
deduce the square-billiard distribution (13). 

In Fig. 2 the wall-collision distribution, which is predicted by Eq.(13), is compared with 
that simulated in square billiard at observation time tobs = IOOTco (for experimental details, 
see Ref. [6]). Analysis for whole-scale temporal evolution was also performed [13]. We there- 



fore infer that the simulated root-mean-square deviation for collision numbers y^A^n^^^i) 
agrees with the stationary prediction given in Eq.(12), starting with times tgbs ^ 30rco (see 
the inset in Fig. 2). 



C. Random Walks in Lorentz Gas Lattice 

1. Infinite Horizon Geometry 

In the case of the scatterer radii < R < L/2, besides a free motion along the open 
corridors in the phase space of the chaotic SB, particles are dispersed by disks and in this 
way are involved in diffusive motion. Let us describe an evolution of a given trajectory, 
moving in the equivalent square-lattice LG, by random walking (see inset A in Fig. 1). 

In a time t, a walker scattered by disks indicates random steps: n{t) — s'^ + s~ + Sy +Sy. 
These steps are given by interceptions of a trajectory with the LG unit cells, that corresponds 
to reflections from the SB walls. More precisely, the walker under consideration does s+ steps 
to the right, s~ steps to the left in the x-direction as well as Sy steps to the down, and Sy 
steps to the up in the y-direction. The resultant particle displacement is therefore given by 

Ar{t) = (4 - s-)i,e, + (4 - s-)iyey. (14) 
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Here and iy are projections of the one-step, free-motion displacement that occurs between 
the two consequent intersections. For random walking in the x-direction, < s+ — s'^ >c= 
within the same approximation, and thus < Ar >c= 0. For the mean squared displacement 
< ArAr >c, one has 

< AV >e= 2il {< s\ >c-< SiSk >c (1 - Sik)] , (15) 

with the help of Eq.(14), where Sik is the Kronecker symbol. Here the mean < il+i^ >c— 2 < 
£l >c— is estimated in the isotropic approximation. The introduced indices i,k — 1, 2, 3, 4 
count distinct random steps: Si = , S2 = s~ , Ss = Sy , and = s~ , which are dynamically 
equivalent: < Sj >c—< Sk >c, for i ^ k. This observation permits one to describe the 

4 

random- wall collisions through the random- walk steps, on the basis of the relation n— X] 

1=1 

This results in the mean Uc and the variance A^ng, namely 



n, 



^{si)c = 4(si)c, AV = 5^(A'si)e + 2 ^(AsiAsfe)c, (16) 

1=1 1=1 i>k 

introduced by the random step fluctuations Asj = Sj — (sj), discussed above in Eqs. (3) and 
(6). 

In the stochastic approximation, the step fluctuations are independent and the last term 
in Eq.(16) is null. Moreover, one can see that the term (A^Si)c = {< >c — < Si >l) is the 
same as in the square brackets in Eq.(15). This leads to Eq.(6), that in combination with 
Eq.(7) provides the desired relation for the variance for collision deviations, namely 



where the characteristic billiard coUision time Tc{R) is given in Eq.(2). 

On the basis of Eq.(17), wc have elaborated a numerical statistical analysis for the wall 
collisions in SB with < R < L/2. As seen from the upper inset in Fig. 3, the standard 
deviation indicates two superdiffusion motion regimes, which are well distinguished through 
the dynamic exponent z{R) < 2. For small scatterers, the /^-independent and, therefore, 
universal regime is manifested by zi — 1.50 ±0.05 (shown by closed squares in Fig.3). Above 
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a crossover radius w 0.35L, the diffusion exponent continuously increase with it! from 1.5 to 
2 (shown by open squares in Fig. 3). A crossover from the universal to a transient dynamics 
starts at Ri = \/2L/4 « 0.35L, the point where a rearrangement of the LG lattice occurs 
with closing of the diagonal Bleher corridors (shown in inset B in Fig. 1). 

We focus on the late-time asymptotic collision behavior. The observed distributions 
have been therefore reestimated in the reduced coordinates proposed in Ref. [3]. These 
coordinates can be formally introduced here by the relations 

D(n) = y/^^D[n(n, t),t] and n(n, t) = ^^^f^, (18) 

^yA^nc{t) 

where n{n,t) stands for the inverse function to n{n,t). The distribution D{n,t), the mean 
nc{t), and the variance A^nc(i) are given in Eqs.(4), (3) and (17), respectively. Physically, 
the reduced collision number n{n,t), as well as its distribution D(n), is expected to expose 
a stationary behavior when the late-time observation condition tobs ^ TciR) is satisfied 
(see also Eq.(7)). As follows from our temporal analysis shown in the insets in Fig. 3, the 
superdiffusion dynamic regime becomes steady starting from the observation times t^^^f^ <; 
50tc. The same can be referred to the steady collision distributions exemplified in Fig. 4. 

As seen from Fig. 4, both the kinds of collision distributions are similar. With growth 
of the scatterer radius, they change their form from characteristic for the ordered motion 
(illustrated in Fig.2 ) to that which tends to the disordered Gaussian motion, shown by solid 
lines in Fig. 4. Furthermore, according to simulation studies carried out within the domain 
0.05 < R/L < 0.35, the late-time distribution functions fall down into the coinciding curves 
Di{n), attributed to the universal superdiffusion regime introduced above by the exponent 
zi — 1.50 ± 0.05. This regime is ensured by the open diagonal and nondiagonal corridors, 
that makes plausible to adopt that the mean collision time is caused mostly by the wall 
reflections. This implies that Td w t^i^ < t^I^ and thus n^^^ > n^^i estimated for R < Ri 
(see Eq.(8)). Consequently, the mean coUision number rid —< ?T'i(x, t) >c~ t/rd^ , defined 
in Eq.(3), can be specified through random walks, namely 
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ni(x,i) = y — ^ , (19) 

j f-cl I'^V 'cl 

with r = i = -y/x^ + defined by the position set x = (x,y,9). Eq.(19) describes the 
scatterer-scatterer collision, which occurs in the periodic LG during time t between two 
scatterers connected by the vector x, with ^^Xj = x; t'^\6i) is defined by Eq.(lO) for a 
given linear trajectory i, which corresponds to the two consequent wall-to-wall collisions. 

On the other hand, the superdiffusion regime under discussion can be treated through 
the two-dimensional Levy jumps, which occur between two scatterers of distance r. There- 
fore, the universal regime can be additionally characterized by the waiting-time probability 
distribution function 

*(r, t) = As{r)5{r - t) with As{r) oc ri~^ (20) 

with 1 < z < 2 and the jump- length distribution As{r) function presented in the long-tail- 
distance asymptotic form. Eq.(20) represent Eq.(39) in Ref. [14], deduced with the help 
of Eq.(39) in Ref. [14] juxtaposed with Eq.(7). In turns, the asymptotic scatterer- collision 
distribution function 

D['\n) = A,[r(n)]^ a ri"^ (nri"^) , with r(n) = nrif^ » 1, (21) 

follows from Eq.(20) and the last relation in Eq.(19). Numerical analysis of Eq.(21) is given in 
the reduced semi- log coordinates in the left inset in Fig. 4. The best fitting for the proposed 
D[^\n) with simulation data results in the derived dynamical exponent zi = 1.5, which is 
the same obtained in Fig. 3 by the variance-number analysis. This justifies the mechanism 
of long-range Levy jumps in the Bleher principal corridors. The universal superdiffusion 
motion is also exposed by numerical analysis for the wall- collision distribution D^\n), 
given in the right plot in Fig. 4. 

The velocity dispersive parameter a{R) was introduced in Ref. [5] to describe a disorder- 
to-order crossover in the weakly open SB in terms of the normal-to- wall- velocity (v_l) pseudo- 
Gaussian distribution function ga{v), with v — v_l/v. The wall-collision statistics could 
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be also introduced [5] through the velocity-dependent random number given by n{v, t) — 
2vt/Tc^\ with the mean 



Jo 



ricait) = / n{v,t)g„{v)dv = and g„{v) = ^^^^ /n^li^'^^ 

(7v27r err(l/2v2(7) 



Here erf(a;) is the standard error function. Two theoretical predictions were proposed for 
the dispersive parameter, namely 



aiR) = '- — , for R<R2 = -, 23 

obtained in the simple "quasi-chaotic" (cr -C 1) approximation [5] and [13] 



a{R) = (7(0)' 



• / 1-2_R/L 

arcsm — 



^c.x.ox..v , ; 1-2R/L 



which takes into consideration the long-living bouncing ball orbits and the square-billiard 
data [13] a{0) = 0.29. In the right inset in Fig. 4, the simulation data for the weakly open SB 
is compared with the theoretical predictions made for a{R). Remarkably, that similarly to 
zi{R), the late-time dispersion parameter is i?-independent, within the experimental error, 
and can be therefore characterized by (Ji{R) = 0.14 ±0.02 for the universal diffusion regime. 



2. Finite Horizon Geometry 



In the case of LG with R > L/2, when all the Bleher corridors are closed (see inset B 
in Fig.l), an application of the central limit theorem for random- walk displacements r[t) 
proves [1] its Gaussian distribution. By taking into account the established in Eq.(6) rela- 
tion between the random displacements and collision numbers, one may expect the normal 
collision distribution, namely 



n,t) 



1 



v/27rA2ne(t) 



exp 



{n-nc{t)y 



, for i? > L/2 



(25) 



2A^n,{t) 

where Uc and A^Uc are the standard mean and deviation of the wall-collision number. 
Experimental justification of the asymptotic Gaussian distributions D^\n) — D2^(n) — 
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(l/27r) exp(— for both the wall and the scatterer coUisions are given in Fig. 5. In 
general, our data on D^^ (n) are consistent with the first observation of the normal scatterer- 
collision distribution reported in Ref. [3]. Furthermore, our short-time analysis provides 
evidence that the Gaussian distribution, associated with chaotic dynamics, becomes steady 
at times tobs ^ 50tc. 

It is noteworthy that no true Gaussian distribution was achieved for the late-time dif- 
fusion coefficient [2], when R tends to L/2 from below. With the aim to clarify a decay 
process of spatial correlations within this regime, we have analyzed the reduced fourth-order 
moments (4 defined as 

^^(^) = ^'^^ - » 1' (26) 

and obtained with the help of Eq.(5). The true Gaussian distribution prescripts ^4 = 3. 
Our temporal analysis given for the case ol R/L — 0.6, results in C,4{t) = 3 ± 0.01 for 
the observation times t > 70tc . The distinct cases of the wall and scatter collisions are 
described by the characteristic times Tc given in Eq.(2). An approximation towards the 
Gaussian distributions, with the increasing of the disorder with R/L, is exposed in the insets 
in Fig. 5. One can see that a chaotic motion in the SB with i? > L/2 is estabhshed by both 
the normal collisions and by the normal diffusion. These universal regimes are observed 
with a good precision, which is guaranteed by, respectively, the higher-order central moment 
C4 = 3.00 ± 0.04 and the diffusion exponent Z2 = 2.0 ± 0.1. 

III. CONCLUSION 

We have discussed the superdiffusive behavior of SB, as well as of the dynamically equiv- 
alent periodic LG, in view of the interplay between the linear and the circular boundary 
geometries. The late-time correlations, driven by elastic reflections from the square walls 
and the dispersed disk, are shown to be steady for a given geometry at observation times 
30-50tc(-R). This justifies a consideration of the statistical distributions for the random wall 
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and scatterer collisions. Besides the higher-order correlation effects, these collisions are char- 
acterized by their mean collision numbers ni"'^ and ■n!i\ respectively. The relative collision 
numbers n'^^ /ni^^ = Tc^^ /t^^ = 2L/t:R follow from Eqs.(2) and (8). 

When the scatterer is absent, the deterministic description results in the ballistic motion 
given by the diffusion exponent Zq — 1. The wall-collision distribution is fiat and asym- 
metric due to a sharp contribution from the bouncing-ball orbits [8] (see Fig. 2). At very 
small scatterer radius, when the wall collisions predominate {nh"^ ^ "n^c^), the late-time dis- 
tribution tends to a smooth pseudo-Gaussian form, asymmetrically shifted from the center 
by the same long-living bouncing-ball orbits (see Fig. 4). When the radius of scatterers 
is relatively small, R < Ri — L-\/2/4, all the principal free- motion corridors remain open 
in the corresponding LG, and the bouncing-ball trajectories evolve freely within the diago- 
nal and nondiagonal corridors. This evolution, described by random walks driven by rare 
collisions with the scatterers, is revealed through the long-distance Levy jumps along the 
diagonal corridors. This stochastic motion is shown to be radius-independent and charac- 
terized by the universal diffusion exponent zi = 1.50 ± 0.05, as well as by the constant 
pseudo-Gaussian disorder-order parameter ai — 0.14 ± 0.02. It seems plausible to associate 
the underlying mechanism of the universal superdiffusion dynamics with the critical behavior 
of spatial correlations [15], with Zcr = 3/2, described by the unbounded orbits moving in the 
diagonal Bleher corridors. When the diagonal corridors become closed, one observes a tran- 
sient motion regime for Ri < R < R2 and characterized by the nonuniversal diffusion with 
Zi < z{R) < 2. For the geometry with R > R2 = L/2 , the scatterer overlaps the boundary 
walls and all the free-motion corridors close. In this case, the disordered motion is due to the 
bounded-orbit normal diffusion with Z2 — 2 and to the normal collisions with n'^^ « n'f^ 
The bouncing-ball motion disappear and the central symmetry of the collision distribution 
is recovered. Despite of the fact that the normal correlations were earlier expected for this 
billiard geometry from rigorous theory [1] and observed in numerical experiments [3,8], we 
have proven that both the observed wall and scatterer collisions are true Gaussian. 

Financial support of CNPq is acknowledged. 
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FIGURES 



Fig. 1. Sinai billiard and the corresponding periodic Lorentz Gas model for different 
boundary geometries. Case A: reduction of a trajectory in the LG lattice to the corre- 
sponding orbit in SB is shown through the billiard wall-to-wall and wall-to-disk collision 
points 1, 2, 3, 5 and 4, 6, respectively. Case B: infinite horizon geometry for < < L\/2/4. 
Examples of the principal free-motion corridors by Bleher [2] : the nondiagonal (a) and di- 
agonal (b) corridors. Case C. finite horizon geometry, R > L/2. Case D: reduction of the 
unbounded trajectory of launching angle 9 in LG without scatterers to the corresponding 
orbit in the square billiard (R — O). 

Fig. 2. Analyses of the wall-collision dynamics in square biUiard. Points correspond 
to simulation data on the collision function observed at to6s — IOOtco and on the temporal 
evolution for the root-mean-square collision- number deviation, with Tco = 7r/4. Solid lines 
are the corresponding theoretical predictions given by Eq.(13) and Eq.(12). 

Fig. 3. Diffusion dynamic exponent against the reduced radius of the scatterer disk in 
Sinai billiard. Points are simulation data derived from the observed billiard-wall collisions 
through the standard deviation (17) and the estimated characteristic time t^\r) given in 
Eq.(2). Inserts: temporal evolution of the standard wall-collision- number deviation for dif- 
ferent infinite {R < L/2, upper inset) and finite ( R > L/2, lower inset) horizon geometries, 
in the log-log coordinates. The solid lines are the best linear fitting of the simulation data. 

Fig. 4. Collision distributions against collision numbers in Sinai Billiard with the infinite 
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horizon geometry (it! < L/2). Points are simulation data for the scatterer (left plot) and 
the wall (right plot) collisions represented in the reduced coordinates defined in Eq.(18). 
The solid lines are the Gaussian distributions D2{n) = (l/27r) exp(— n^/2). Left inset: 
analysis of the late-time scatter-coUision function for long-jump-scatterer collisions. Points 
are simulation data for the case R — 0.25L and tabs — '200tc^\ Sohd hne is the best fitting 
with the prediction given in Eq.(21), with D^\n) = 0.5|n|~^^/^. Right inset: velocity- 
dispersive parameter against reduced scatterer radius. Points are simulation asymptotic 
data [13] for the weakly open SB; dashed and sohd lines are theoretical predictions given in 
Eqs. (23) and (24), respectively. 

Fig. 5. Observation of the Gaussian collision statistics in Sinai Billiard with R > L/2. 
Points are simulation data for the scatterer (left plot) and the wall (right plot) coUisions 
observed at tabs — 200tc for distinct geometries indicated in the legend. The lines are the 
same as in Fig. 4. In the left and right inserts, the analyses of the fourth moment (defined 
in Eq.(26)) for the scatterer and wall coUisions are given at tabs — 200tc, respectively. Points 
are simulation data and the dashed line corresponds to (^4 = 3. 
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